Sintering behaviour of two roughened crystals just after contact 
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We consider two spherical, roughened crystals with approximately isotropic surface free energy 
which are brought into contact and begin to sinter. We argue that the geometry immediately post- 
contact is two dimensional and Cartesian and can be approximated by the evolution of a slot-shaped 
cavity. On this basis, we construct travelling wave solutions for the crystal shape in the limits of bulk 
diffusion limited and surface diffusion limited kinetics. These solutions are then used to calculate 
novel scalings for the neck size as a function of time t after contact: We predict that neck size 
is proportional to t^^'^ for the bulk diffusion limited case and (following a single pinch-off event) 
approximately proportional to t^''^ for the surface diffusion limited case. 

PACS numbers: 81.10.Aj 
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I. INTRODUCTION 

A crystal in equilibrium against its vapor or a melt 
eventually reaches a shape which can be obtained by the 
Wulff construction Let the free energy per unit area 
of a plane surface of the crystal which is perpendicular 
to a unit vector fi be given by 7(n). We define the points 
p(fi) oc n7(n) for some fixed constant of proportionality 
which eventually sets the size of the crystal. Next we 
construct a plane through each p(n), perpendicular to 
n. The Wulff shape of the crystal is the inner envelope 
of all such planes. 

From the terrace-ledge-kink (TLK) (3,01 model we ex- 
pect the surface free energy of a crystal to have cusp-like 
minima near to the crystallographic symmetry directions: 
For a plane oriented at a small angle 9v to a crystallo- 
graphic axis (i.e. a 'vicinal plane') and at absolute tem- 
perature T, the surface free energy per unit area behaves 
like [1 



7(0.,T)«7o(r) + 



a(T) 
dl 



tan( 



(1) 



where 70 is the free energy per unit area of a molecularly 
flat surface with 9v = 0, while a is the ledge energy per 
molecule and dm is a molecular distance. 

Using the Wulff construction, this leads to the emer- 
gence of crystal facets at equilibrium, in which a macro- 
scopic portion of the crystal surface is molecularly flat, 
save for isolated single molecule islands and surface va- 
cancies 0. 

However, at a particular temperature (the 'equilibrium 
roughening' or 'surface melting' temperature), the free 
energy barrier to formation of new molecular islands on 
the facet vanishes; the facets become rough on a molec- 
ular scale, so that the relevant cusp in the surface free 



energy disappears along with the facet, to be replaced 
by a macroscopically smooth, rounded crystal surface 
[3, El- This roughening transition is different for dif- 
ferent crystallographic symmetry directions; for example 
three different temperatures have been observed for He- 
lium crystals (TJ , while for ice, the basal facet persists up 
until melting, while prism planes roughen at — 2°C when 
against vapor Q or — 16°C when against water under 
pressure [3| . If the crystal surface is not at equilibrium, 
but instead growing, this can move the roughening tran- 
sition to lower temperatures [l3, Ell • 

In this investigation, we start from two assumptions: 
First, we assume that all the relevant crystallographic 
planes are roughened (there may however be other direc- 
tions not involved in the analysis described below, which 
are faceted). Second, we make the approximation that 
the roughened portion of a crystal surface has a fairly 
isotropic surface free energy per unit area 7 (a typical 
approximation for phase-field studies of dendritic growth 
or solvability theory [ul). 

When a roughened crystal does not have its Wulff 
shape, then the surface free energy will not be equal ev- 
erywhere. The chemical potential /i of a molecule at a 
curved surface differs from the value fiQ it would have for 
a flat surface by the well known Gibbs- Thomson equa- 
tion, which in linearized form gives 



(2) 
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Here fl^ is the molecular volume and k is the mean cur- 
vature of the surface (positive if the surface is convex) . 

Molecules will therefore have a tendency to leave con- 
vex portions of the surface and attach to concave (or 
less convex) regions. In doing so, they may either travel 
through the bulk of the material outside the crystal, 
in a process of evaporation-condensation (or solution- 
precipitation for a crystal against a melt), or they may 
diffuse along the crystal surface [13]. 

This leads to the crystal surface having a normal veloc- 
ity Vm which will in general be the sum of a component 
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v!^ from molecules arriving by diffusion through the bulk 
and from those arriving along the surface. In sections 
ITTl and mil we consider limiting cases where only one of 
these mechanisms predominates. For vapor phase sin- 
tering, it is likely that surface diffusion will dominate, 
while for liquid phase sintering bulk diffusion may be the 
important mode IJ]. The case where crystal facets in- 
terfere [l^ has been the subject of more recent work, but 
is not our concern in this paper. 

Consider now two crystals which are close to their equi- 
librium shapes (and therefore close to being spherical 
at least in the region of interest). If these crystals are 
brought into gentle contact, the curvature near the con- 
tact point will be large and negative, and therefore a 
neck will appear and rapidly grow. There will usually be 
a grain boundary coincident with this neck, but in con- 
trast to Ref. |14| . where equilibrium is approached from 
an initial dihedral angle of tt radians, in the case of two 
sintering crystals, the initial dihedral angle is 0, and we 
no longer expect grain boundary grooving [l^ . 

The initial stages of neck growth occur rapidly, but 
nevertheless can be of importance in the rheology of crys- 
tal slurries, where contacts may be formed and broken in 
rapid succession. Relevant examples might include the 
flow of lava [l^ or of brash ice [l7| . 

The initial growth of necks between two crystals has 
been analysed previously, for example in Ref. [l^ where 
a i^^^ time dependence for neck radius as a function of 
time t after contact was found for the case of bulk diffu- 
sion limited growth. 

For the development of this paper, we note that there 
is a geometrical simplification which occurs in the earli- 
est stages of neck growth and which appears to have been 
overlooked in the literature: If we consider two crystals 
just after they have touched (Fig. [T]) the parts of the 
surface near the contact point (except for the very high 
curvature neck) are very nearly parallel to one another. 
This suggests that an interesting geometry for investi- 
gating the early stages of neck growth is a 2 dimensional 
slot-shaped cavity, illustrated in Fig. [T](c) . The closer in 
time we are to first contact the better this is as an approx- 
imation. Furthermore, if the radius of the neck is large 
compared to the separation of these two nearly parallel 
surfaces (as is indeed the case at early times), then the 
difference between axi-symmetric cylindrical coordinates 
and plane Cartesian coordinates becomes negligible. 

In sections [Tl] and IIIII we analyse a slot-shaped cavity 
of this kind, and then compensate for the remaining dif- 
ferences between cylindrical and plane Cartesian geome- 
try during structure evolution by invoking conservation 
of volume. In this way, we obtain the geometry of the 
crystals in the earliest moments after contact, and also 
predictions for novel scaling behaviour of the neck size as 
a function of time after contact. 
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FIG. 1: Image (a) shows schematically two spherical crystals 
that are in contact and have started to sinter together, form- 
ing a narrow neck of radius r„ — 0.02 (where all distances 
are scaled by a suitable capillary length). The crystals are 
axi-symmetric about an axis up the page. In (b), a magnified 
portion of the neck region is depicted, illustrating that at the 
very earliest stages of sintering, the surfaces of the two crys- 
tals are nearly parallel close to the neck. Image (c) shows a 2 
dimensional slot-shaped cavity, with exactly parallel surfaces, 
which we use as an approximation to the ajci-symmetric ge- 
ometry in (b). Image (d) shows a 2d slot-shaped cavity of 
finite length, which is more appropriate for simulations. 



II. BULK DIFFUSION LIMITED CASE 

Consider the case where growth is limited by bulk dif- 
fusion to and from the crystal surfaces (attachment ki- 
netics being unimportant because the crystals are rough- 
ened) . 

We assume that the mean free path for bulk diffusion 
is small compared to the feature sizes of interest. If 0(r) 
is the volume fraction in the vapor or melt of the type of 
molecules that can make up the crystal, then this quan- 
tity obeys a diffusion equation with a diffusivity D. 

We make the further approximation that the transients 
of this equation have died out, so that the concentration 
field (j> obeys Laplace's equation 



(3) 



This requires that the surface normal velocity w„ should 
be small compared to D divided by a typical feature size. 

The normal velocity w„ of the crystal surface (which we 
denote by the set of points S) is obtained by conservation 
of volume for the crystallizing material: 



(4) 



rGS 



where n is an outward unit normal vector to S. 

We note in passing that we can apply the divergence 
theorem directly to Eqs. (jH) and with the result that 
the volume of an enclosed cavity is conserved under this 
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FIG. 2: Numerical solution for the bulk diffusion limited evo- 
lution of a parallel-sided 2 dimensional slot-shaped cavity in a 
crystal. The initial surface is shown in bold, and subsequent 
points in time show the formation of dumbbell ends to the 
curve. Although not shown, the final equilibrium shape must 
be a circle at the origin containing the same area as the initial 
curve. 



evolution, while for a discrete crystal in an open volume 
of vapour or melt, one must also consider the concentra- 
tion at infinity before it is possible to say anything about 
conservation of crystal volume. 

Let (j>o{T) be the concentration of molecules in the melt 
or vapour which is present at equilibrium against a flat 
crystal surface at temperature T, then perturbing this 
linearly using Eq. ([5]) we find that the actual concentra- 
tion at the surface of a curved crystal is given by 



flyTjK d(j)Q 

L dT ' 



(5) 



where L is the latent heat per unit volume of fusion or 
evaporation. 

Eqs. pi4|5|) determine the evolution of the system in an 
exactly analogous manner to the LSW theory of Ostwald 
ripening jH^. 

Suppose the co-ordinates of a point in space are de- 
noted by (x,y,z) and time after contact by t. We 
define a capillary length Ic = jfly/L and dimension- 
less space, time and concentration variables x = x/lc, 
t = iTD{d4o/dT)/ll and V' = (0 - (f>o)/{Td(Po/dT). 

Eqs. pi4l5p then take the following simple form in 2 
dimensional Cartesian co-ordinates: 



\ dx^ dy"^ 



ip [x,yo{x,t),t] = 



d'^ya/dx^ 

1 



dm 
dt 



dyo di^ 
dx dx 



(dyo/dx) 
dy 



y=va 



(6) 
(7) 

(8) 



V=Vo 



where y = yo{x) is the crystal surface S. 

If we solve Eqs. (|6I7I8|) numerically, using a simple ex- 
plicit timestepping scheme for an initial slot-shaped cav- 
ity, we see from Fig. [2] that dumbbell-shaped ends start 
to appear as the cavity evolves towards its eventual equi- 
librium shape, which must be a circle. 

This behaviour suggests that there may be a travelling 
wave solution to Eqs. (|6l7l8p which propagates at con- 
stant velocity v and without change of shape. This so- 
lution would consist of a teardrop-shaped cavity formed 



from the evolution of a semi-infinite slot-shaped cavity 
in which the parallel surfaces almost, but do not quite, 
touch. To picture this slot, imagine the initial cavity of 
Fig. [21 but with two changes: First, let it extend from 
2: = — ootoa:; = 0, rather than the range x = —56 to 
a; = 56 (as shown in Fig. [2]). Second, let it be very nar- 
row, rather than extending between y — —6 and y — 6 
(as shown in Fig. [5]). 

If such a solution exists, then it is natural to tackle 
Eqs. (|6l7l8p through complex analysis: Let 



(9) 



then if ^ is an analytic function of z (with an addi- 
tional time dependence), Eq. © is satisfied as an iden- 
tity, which simplifies the problem greatly. 

By balancing leading powers, we find that ip{z, t) must 
have a branch point and the solution must be (up to 
translational symmetry) of the form 



V'(z,i) EE w^/^V) v'^/'^{z + vt);a 
yo{x,t) = v'^/^yo v'^^'^{x + vt);a , 



(10) 



(11) 



where a is a real parameter which represents a constant of 
integration in Eqs. (|6l7l8p . Physically, a determines the 
curvature of the surface near to the tip of the teardrop- 
shaped cavity (the origin in Fig. [3]) . 

After a considerable amount of back-substitution, one 
obtains a series solution for the analytic function ^{z; a) 
and the real function yoix; a): 



7/ \ -1/2 °' 1/2 16 5 

w(z\a) = az ' z ' z ; 

' ^ 18 15 1224 
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The series solution for j/o {x) does not converge quickly 
enough to be useful. However, Eq. for ip{z) does 
converge quickly in the region of interest, and using this 
solution, the profile yo(x) can be drawn numerically, us- 
ing Eq. ([71). The results for two different values of the 
a parameter are shown in Fig. [31 and a graph of the en- 
closed area in Fig. [4l 

Although the travelling wave solution just presented is 
an idealized case, an approximate analytical solution of 
this kind can also be used to construct the rate of neck 
growth in the original problem of two sintering spheres. 
This is done in the following manner: 

From Eqs. (jlOlllI) we see that under magnification the 
velocity of the travelling wave is inversely proportional 
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FIG. 3: Profiles of the traveUing wave with v — 1, for (a) 
a = 0.468703 and (b) a = 0.3. The cavity moves to the left 
with speed v = 1 without changing shape. The thin (in this 
limit, zero thickness) parallel sided slot is shown as the line 
from a; = along the negative real axis (which is the branch 
cut of Eq. ([12|. When a ^ 0.468703, the profile has a cusp 
on the right hand side. This represents the (more realistic) 
solution when the thin slot terminates at a grain boundary 
between two crystals, rather than being a cavity in a single 
crystal. 



l7t 



0.75jt 



0.5jt 



0.25;t 



Oji, 



dihedral angle at grain boundary 
area within cavity 



- 1.5 



1 O 



-0.5 



0.1 



0.2 



0.3 



0.4 



a 



FIG. 4: Left hand scale and solid curve shows a plot of the 
dihedral angle 9o at the grain boundary (if present) as a func- 
tion of the parameter a. The right hand scale and dashed 
curve shows the enclosed area of the travelling wave cavity, 
defined by I{9o) = / 2yo [x;a{6o)]dx. Both plots are for a 
(non-dimensional) wave velocity v = 1. 



to its enclosed area, and is independent of the separation 
of the parallel surfaces in the initial slot-shaped cavity 
(provided this separation is very small). However, if the 
separation between these surfaces is not zero, then the 
teardrop shaped cavity illustrated in Fig. [3] must grow in 
size as it moves, in order to maintain conservation of en- 
closed area (or volume in the three dimensional case), as 




FIG. 5: Schematic illustration of where the travelling wave 
solution of Fig. [3] is expected to occur during the sintering of 
two spherical particles. Left hand illustration shows the two 
crystals soon after contact. Right hand panel shows a mag- 
nified portion of the neck region, with the teardrop travelling 
waves forming the edges of the neck. 



discussed above. We note that for the case of a neck grow- 
ing between two touching spherical crystals, the teardrop 
cavity which will form (shown schematically in Fig. [5]) is 
not a completely enclosed volume, but will be nearly so 
at early times. 

Now suppose we have two spherical crystals with ini- 
tial (non-dimensionalized) diameters di and ^2, which 
are brought into contact. Let the radius of the neck at 
(scaled) time t after contact be r„(i). The volume of ma- 
terial [to leading order in r„(i)/(ii] that has been added 
to the crystals in order to make this solid neck is given 
by 



^(rfi+d2)r^(t)/(2did2). 



(14) 



If we assume that the neck region forms a nearly com- 
pletely enclosed volume, then the volume of crystal in this 
region must be conserved (using the argument from the 
divergence theorem above). The material used to make 
the solid neck must therefore come from the toroidal cav- 
ity formed by the travelling wave (illustrated schemati- 
cally in Fig. [S]). 

We can estimate this volume using Eq. (ITlT) and Pap- 
pus' second centroid theorem 22] (that the volume gen- 
erated by rotating a plane figure around an axis is the 
product of its area and the distance moved by its cen- 
troid). The result (again to leading order) is 



2Trr^{t)v'^I{eo), 



where 



I{8o) = / 2yo{x;a)dx 



(15) 



(16) 



is the enclosed area, which is a function of the equilibrium 
dihedral angle Oq at the grain boundary. The dihedral 
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angle in turn is set by the grain boundary surface energy 
per unit area 7g b through [ij] 



27Cos(6'o/2) = 7GB. 



(17) 



By equating the two volumes of Eqs. p4)) and (IT5|) . we 
obtain the final result for the growth velocity of the neck 
at small times: 



so 



dt 

that r„(t) cx t^/'^. 



ri di + d2 



(18) 



III. SURFACE DIFFUSION LIMITED CASE 

If recrystallization is limited by the rate of surface dif- 
fusion instead of bulk diffusion, then (still for the rough- 
ened, nearly isotropic 7 case), the normal growth velocity 
of the crystal surface is given by [T^ 



(19) 



where Ug is the number of molecules per unit area at the 
surface available to diffuse, Dg is a surface diffusivity, ks 
is Boltzmann's constant and V| is the surface Laplacian 
(also known as the Laplace-Beltrami operator [23j). 

Under the action of Eq. the volume of the crystal 
is exactly conserved, which follows from the equivalent 
of the divergence theorem for V| 23| . 

In order to proceed, we non-dimensionalize in a sim- 
ilar manner to Sec. |TT1 Let the co-ordinates of a point 
in space be denoted by {x,y,z) and time after contact 
by i. We define a capillary length Ig = jilv/ksT and 
dimensionless space and time variables x = x/lg and 
t = iflyUsDs/l^. 

The equation for surface evolution then takes the non- 
dimensionalized form 



(20) 



If we again perform a simple, explicit numerical sim- 
ulation for the evolution of a slot-shaped cavity in 2 di- 
mensions, we obtain the results shown in Fig. [6l 

This problem has in fact been studied before [13, HI] , 
in the context of the 'inverse' case, where the authors 
study the evolution and breakup of a thin plate. How- 
ever, since the evolution is determined by the surface 
only, this does not affect the equations. The formation 
of dumbbell shaped ends in Fig. [5] is then the analogue of 
the first stages of edge spheroidization of a plate (2J, [2^ . 

Once more, this figure suggests that it is worth seek- 
ing a travelling wave solution to the evolution equation, 
which moves at constant velocity v without change of 
shape. To do this, we note that if Q is the angle that the 
curve representing the crystal surface makes with the y- 
axis, s is the dimensionless distance along this curve and 
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FIG. 6: Numerical solution for the surface diffusion limited 
evolution of a thin, parallel-sided slot-shaped cavity in a crys- 
tal. The initial surface is shown in bold, and subsequent mo- 
ments in time show the formation of dumbbell-shaped ends 
to the curve, before the final equilibrium shape is reached, 
which is a circle enclosing the same area as the original curve. 
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FIG. 7: Travelling wave solution for the equation of surface 
diffusion limited evolution [Eq. (1201) ]. The entire curve moves 
to the left at speed v — 1. This is a purely mathematical 
solution to the problem, as the curve intersects itself and so 
cannot represent a real crystal boundary. The total (algebraic 
and non-dimensionalized) area enclosed by the curve is 2.536 



V is the dimensionless speed of propagation, then the rel- 
evant equation to solve is [from Eq. (|20|) ] 



V cos 6 — 



d^ejs) 

ds3 



(21) 



with boundary condition 6(0) = rf^9(0)/ds^ — 0. 

Fig. [7] shows a numerical solution corresponding to 
V = 1, while the form of Eq. (pij) shows that under 
magnification, the velocity v is proportional to the en- 
closed area to the power —3/2. Although a formal solu- 
tion to Eq. (PT|) , the curve is self-intersecting; indeed the 
assymptotic form for large negative x is 



yoix) ±e 



x/2 



COS 



{xV3/2). 



(22) 



This self-intersecting shape strongly suggests that un- 
der the evolution equation [Eq. ([20|) ]. a long slot-shaped 
cavity will repeatedly pinch off to form a string of equally 
spaced circular cavities in its wake. This behaviour is 
indeed seen in the numerical simulation of a very long 
cavity in Fig. [8l 

From the simulation results of Fig |S] and the scaling 
properties of Eq. (PT|) , we see that for a long slot-shaped 
cavity with a dimensionless separation w between the 
opposite, parallel crystal surfaces, the average (dimen- 
sionless) velocity of the end of the slot will be 



(v) « 0.0263iy 



-3 



(23) 



and the radius of each circle left behind in the wake is 
approximately 5.211;;. 
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FIG. 8: Snapshots in time (displayed with vertical offset for 
clarity) for the evolution of an initial long slot-shaped cav- 
ity under surface diffusion limited evolution [Eq. (|20|) ]. The 
initial cavity is shown in (a) and has a dimensionless width 
w = IB/tt. Configuration (b) is the first pinch-off event, at 
dimensionless time t — 2.17x 10®. Image (c) shows the second 
pinch-off event at t — 4.34 x 10® and (d) is the final, equilib- 
rium configuration. Curves are drawn only for a; > 0, and are 
in fact symmetric under reflection about x — 0. 



For a pair of spherical crystals just after contact and 
evolving according to Eq. ((20)) . it would be tempting to 
use Eq. (|23|) . replacing w virith the separation of the undis- 
turbed spherical surfaces. This would give an average 
rate of growth of the radius r„ of the neck proportional 
to r~^. However, this would be incorrect, as conservation 
of crystal volume is not taken into account properly. 

Instead, for the axi-symmetric case of two spherical 
crystals each of (dimensionless) diameter d, we impose 
conservation of crystal volume in the following manner: 
First, assume that pinching-off produces a series of con- 
centric tori with major radii given by the set {ri} with 
« S {1, 2 . . .}, and minor radii by the set {pi} given by 



5.21w{n) EE 5.21 d - Jd^ - 4r 



(24) 



{w{r) being the separation of the surfaces of the origi- 
nal spheres at distance r radially from the initial contact 
point). 

Conservation of enclosed volume then allows us to 
write down the distance from the axis of symmetry at 
which the i'th torus lies 1221: 



2'Krw{r)dr « 2'Kri'K [5.2lK;(ri)]^ . 

However, Eq. (^5]) has only one solution, namely 
d 



877(5.21)2 



0.00147d. 



(25) 



(26) 



We therefore expect the growing neck to leave behind one 
tiny toroidal cavity, before proceeding to grow without 
pinch-off events. We conjecture that this toroidal cavity 



will subsequently break up into a ring of spheres by a 
similar mechanism to the breakup of a cylinder p6l [ or 
edge spheroidization of a plate [2J|. 

If, after this single pinch-off event, the profile still 
maintains a dumbbell-shaped end, similar to Fig. |S1 and 
containing most of the volume of the vapor or melt in the 
neck region, then (again by conservation of enclosed vol- 
ume and using the algebraic area in the curve of Fig. [7]) , 
the neck radius r„ should grow as 



drn jt) 
dt 



2 X 2.536r„d 
r4 - (0.00147d)^ 



3/2 



(27) 



so that very approximately, r„(i) oc t^^^. 

We note that grain boundaries are also easily incorpo- 
rated into this formalism, by changing the value of 0(0) 
in Eq. ((2T|) to introduce the relevant dihedral angle at the 
origin in Fig. [T] Just as in section [III the grain bound- 
ary runs parallel to the slot (and outside it), while in 
Ref. [iJl the geometry is very different: grain boundaries 
run perpendicular to the plate (and inside it). 



IV. CONCLUSIONS 

For two sintering roughened crystals, the observation 
that the geometry near the contact point in the first mo- 
ments after contact is two dimensional (and very similar 
to a slot-shaped cavity) is a considerable simplification 
for analysing the problem at these early times. 

We are able to bring to bear complex analysis in the 
case of diffusion limited crystal evolution, and simple nu- 
merical approaches for the surface diffusion limited case. 

The result is a novel scaling for the power law for 
growth of the neck in the diffusion limited case [namely 
r„(t) (X t^/'* from Eq. ([T5|) ]. and the prediction of a pinch- 
off event [Eq. (|26p ] followed by approximate power law 
growth [Eq. (P7|) ] for the surface diffusion limited case. 

We also believe that the travelling wave solutions piv- 
otal to these results (namely Figs. El and [7]) have some 
aesthetic appeal. 
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